A PARABOLIC FREE BOUNDARY PROBLEM WITH 
BERNOULLI TYPE CONDITION ON THE FREE BOUNDARY 



JOHN ANDERSSON AND GEORG S. WEISS 

Abstract. Consider the parabolic free boundary problem 

Au -dtu = in {m > 0} , | Vm| = 1 on d{u > 0} . 

For a realistic class of solutions, containing for example aH limits of the singular 
perturbation problem 

Aue — dtUe = f3e{ue) as £ — ♦ 0, 

we prove that one-sided flatness of the free boundary implies regularity. 
In particular, wc show that the topological free boundary d{u > 0} can be 
decomposed into an open regular set (relative to d{u > 0}) which is locally a 
surface with Holder-continuous space normal, and a closed singular set. 
Our result extends the main theorem in the paper by H.W. Alt-L.A. Caffarelli 
(1981) to more general solutions as well as the time-dependent case. Our proof 
uses methods developed in H.W. Alt-L.A. Caffarelli (1981), however we replace 
the core of that paper, which relics on non-positive mean curvature at singular 
points, by an argument based on scaling discrepancies, which promises to be 
applicable to more general free boundary or free discontinuity problems. 



1. Introduction 

The parabolic free boundary problem 

(1.1) Aw ~dtu = Q in {u > 0} , |Vu| = 1 on d{u > 0} 

has originally been derived as singular limit from a model for the propagation 
of equidifFusional premixed flames with high activation energy ([HI); here u = 
X(Tc — T) , Tc is the flame temperature, which is assumed to be constant, T is 
the temperature outside the flame and A is a normalization factor. 
Let us shortly summarize the mathematical results directly relevant in this context, 
beginning with the limit problem 1)1. in the brilliant paper J], H.W. Alt and L.A. 
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Caffarelli proved via minimization of the energy J(|Vup + X{u>o}) ~ here X{u>o} 
denotes the characteristic function of the set {u > 0} - existence of a stationary so- 
lution of in the sense of distributions. They also derived regularity of the free 
boundary d{u > 0} up to a set of vanishing n — 1-dimensional Hausdorff measure. 
By |12j existence of singular minimizers implies the existence of singular minimizing 
cones. L.A. Caffarelli-D. Jerison-C. Kenig showed that singular minimizing cones 
do not exist in dimension 3 (IHI)- Moreover it is known that singular minimizing 
cones exist for ?^ > 7 Non-minimizing singular cones appear already for 7i = 3 

(see example 2.7]). Moreover it is known, that solutions of the Dirichlct problem 
in two space dimensions are not unique (see example 2.6]). 
For the time-dependent Hl.lf) . both "trivial non- uniqueness" (the positive solution 
of the heat equation is always another solution of ULII) ) and "non-trivial unique- 
ness" (see Oni) occur. Even for flawless initial data, classical solutions of 
develop singularities after a finite time span; consider e.g. the example of two 
colliding traveling waves 

^ 2) u{t,x)= X{x+t>i}icxpix + t - 1) - 1) 

+ X{-..+t>i}(exp(-x + t - 1) - 1) for t e [0, 1) 

(see Figure P). 



t = 1 



t = 




X = —1 X = X = 1 



Figure 1 . Colliding traveling waves 



There are several approaches concerning the construction of a solution of the time- 
dependent problem, all of which are based in some form on the convergence of the 
solution Me of the reaction-diffusion equation 

(1.3) Au, - dtu, = f3,{u,) 

to dJ as e ^ 0; here Pe{z) = i/3(|) , G C^i[0, 1]) , /3 > in (0, 1) and / /3 = i . 
L.A. Caffarelli and J.L. Vazquez proved in jTj uniform estimates for (|1.3|) and a 
convergence result: for initial data that are strictly mean concave in the interior 
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of their support, a sequence of e-solutions converges to a solution of in the 
sense of distributions. 

Let us also mention several results on the corresponding two-phase problem, which 
are relevant as solutions of the one-phase problem are automatically solutions of 
the corresponding two-phase problem. In [H] and 0], L.A. Caffarelli, C. Lederman 
and N. Wolanski prove convergence to a barrier solution in the case that the limit 
function satisfies {u = 0}° = . 

Then, there is the convergence to a solution in the sense of domain variations 
which seems to contain more information than the barrier solutions in [S] and 0]. 
For more general two-phase problems see Domain variation solutions play an 
important rule in this paper and will be discussed in more detail in Section |21 

Here let it suffice to say that domain variation solutions are pairs (it, x) where the 
order parameter x shares many properties with the characteristic function X{u>o} 
but does not necessarily coincide with it. By all limits of the singular pertur- 
bation problem H1.3() are domain variation solutions, so all results in the present 
paper hold for all limits of H1.3(l . 

Our main result Theorem 18.41 states - leaving out inessential assumptions - that 
if (0, p"^) is a point on the topological free boundary and if the set {x > 0} is flat 
enough, i.e. 

x(x, t) = when (x, t) e Qp and .t„ > ap, 

for some cr < cto (see FigurelSJ, then the free boundary Qp/4r)d{u > 0} is a surface 
with Holder-continuous space normal. 
As a consequence we obtain that the regular set is open relative to d{u > 0} 



{x,t) = {0, 1) 


















^/x = 








Figure 2. One-sided flatness in the case p = 1 



4 



J. ANDERSSON AND G.S. WEISS 





X > 


X > \ 




N^^^egular^free boundary 1 
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x>0 


X>0 





Figure 3. Example of the set of regular free boundary points (stationary) 
(Corollary iSl cf. Figure 0). 

Note that even in the stationary ease our result extends the result in ^ as our 

assumptions do not exclude degenerate points or cusps close to the origin (excluded 

by the definition of weak solutions ^ 5.1]), our result does that. 

In the proof of our result we use ingenious tools developed in Q]: We prove that 

flatness on the side of {x = 0} implies flatness on the side of {% > 0} which 

in turn yields uniform convergence of an inhomogeneously scaled sequence of free 

boundaries. 

However we replace the core in the method of H.W. Alt-L.A. Caffarelli, relying on 
non-positive mean curvature of d{u > 0} at singularities, by a method based on 
scaling discrepancies (Proposition!^^. This original component gives hope that the 
method may now be applicable to more general free boundary or free discontinuity 
problems, in particular two-phase free boundary problems. 

2. Notation 

Throughout this article R" will be equipped with the Euclidean inner product 
X ■ y and the induced norm |a::| , Br^xo) will denote the open n-dimensional ball of 
center xq, radius r and volume r" ujn , Bl.(0) the open n — 1-dimensional ball of 
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center and radius r, and a the i-th unit vector in R". We define Qrixo,to) :~ 
Br{xo) X {to — r"^ , to + r'^) to be the cyhnder of radius r and height 2r^ , Q" (ccq , to ) : = 
Br{xo) X (to - f^,to) its "negative part" and T!,r(to) := x (to - 4r^,to - r^) 
the horizontal layer from to — 4r^ to to — r^. Let us also introduce the parabolic 
distance pardist((t, x), A) := inff^ ^J\x — yP + |t — s|. Considering a function 
(j) G i?[Q^(R"; R") we denote by div := X]"=i ^i<?^j the space divergence and by 



D(j) 



... dnCpl 



di4>n ... (?„ 



the matrix of the spatial partial derivatives. 

Given a set A C R" , we denote its interior by A° and its characteristic function 
by xa ■ In the text we use the n-dimensional Lebesgue-measure and the m- 
dimensional Hausdorff measure 7i™. When considering a given set A C R", let 

OmA := {x e R" : hmsup — \ )' . > and limsup — \ )' . > 0} 

be the measure-theoretic boundary of A, let d* A := {x E R" : there is ^{x) G 
dBi{0) such that ?'""/b^(^) \xa - X{y:{y-x)-,y{x)<o}\ ^ as r ^ 0} (by [H Corol- 
lary 5.6.8] d*A coincides H"~^-a.e. with the reduced boundary of a set of finite 
perimeter defined in |14l Definition 5.5.1]), and let ly : d*A — dBi{0) denote this 
measure theoretic outward normal to dA. We shall often use abbreviations for in- 
verse images like {u > 0} := {x e O : u{x) > 0} , {a;„ > 0} {x G R" : x„ > 
0} , {.s = t} := {{s,y) G R"+^ ; s = t} etc. as weh as A{t) A n {s = t} for 
a set A C R"+^, and occasionally wc employ the decomposition x = {x',Xn) of a 
vector X G R" as well as the corresponding decompositions of the gradient and the 
Laplace operator, 

Vu — (V'u, dnu) and Au = A'u + dnnU . 



Finally, C^'^ :— H^"'^ denotes the parabolic Holder-space defined in [H]. 

3. Notion of solution and Preliminaries 

In this section we gather some results from QT]. As degenerate points are un- 
avoidable in the parabolic problem (see the introduction of for examples), an 
extension of the weak solutions in pP does not seem to be the right choice. Instead 
we use the solutions of [111 Definition 6.1], which are, roughly speaking, solutions in 
the sense of domain variations. The advantage is that the class of solutions defined 
in ^2 Definition 6.1] is closed under the blow-up process. Moreover, all limits of 
the singular perturbation problem discussed in |2] are domain variation solutions 
and satisfy ^2 Definition 6.1] (see [111 Section 6]). Let us recall the definition of 
solutions and the monotonicity formula used therein: 
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Theorem 3.1 (Monotonicity Formula, cf. ^ Theorem 5.2]). Let {xo,to) e R" 
(0,00) , T;'{to) = R" X {to - Ar^,to - r^) , < p < a < ^ and 

G(.o,to)(a;,0=47r(to-0|47r(to-t)r^"' (" 4(^~ ) ■ 

Then 

satisfies the monotonicity formula 

*(:EO,to)(f^) - *(xo,to)(p) 



X 



> 



T-(to) 



Definition 3.2 (cf. m Definition 6.1]). We cafi (u,x) a solution in fio := R" x 
(0, 00) (in which case we set r 0) or ili := R" x (—00, 00) (in wliicli case we set 
T := 1), if:^ 

i)ue cll{nr)nc^{nrn{u > o})nHlf^inr) and x e LH(-^-R,^);S^^(-Bii(o))) 

for each R £ (0, 00) . For each i? S (0, 00) and 5 G (0, 1) there exists Ci < 00 such 
that for (9r(a;o,io) C n (5i{(0) 



/ |Vx| < Ci 

Qr(a:o,to) 

/" |9twp < Ci r", and 



'B,(xo)x(to+Sir2,to + S2r2) 

for < 5*1 < 5*2 < 00; here the mollifier {4's)5e(o,i) should be non-negative and 
satisfy 4>s{-) ^ 4> e C°'i(R") , / = 1 and supp C Bi(0) . 

Moreover, x G {0, 1} a.e. in and X{ji>o} X s-C- in ■ 

2) The solution u satisfies the monotonicity formula Theorem 13.11 (in the case of 
T = 1 for (.To, to) e R"+^ and a e {Q,oo)). 

3) = /" f [-2dtuVu-^ + (|Vup + x)div^ - 2\/uD^\/u] 

for every ^eC°'^(f7^;R") ■ 

4) The solution u is non-negative. 

5) The solution u attains the initial data w° S Cq'^(R") in Lj^q^(R") in the case 
that T = . 

6) For each k > there is 5 > such that Qrixo,to) C fir and || uixo+rx^+r^t) _ 
^|2;«||lco(Qi(o)) < S imply <1 + k. 
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7) For 6 G (0,1) , e C','\{\y\^ + < S^}) , Ur{y,s) := "(*o+^^^.^o+.y) 
Xriy, s) := + fy, to + the following holds: 

a) f \{\/xr-x + 2tdtXr)*H 

JQp{xi,ti) 



for -S < ti < -T < , S + p'^ < ^ , \xi\ < Z and, in the case of r = , io - 
2r^{-ti+ + S) > 0. 



b) / li^^Xr-O^H < CiS) / |Vu,-ei 

for C e 5Bi(0) , ti < and, in the case of r = , to - r'^{-ti + {p + VS)'^) > . 

dt{{\\i'Ur\'^ +Xr)*(l)s)it,Xo) < I I 2dtUr{t,z)\/Ur{t,z)-\/(j)s{xo- z)dz 

for — oo < ti < t2 < oo and, in the case of t = , to + r^ti > . 

Remark 3.3. As the function x is defined only almost everywhere, all pointwise 
equalities/inequalities involving x should be understood as equalities/inequalities 
that hold almost everywhere with respect to the Lebesgue measure. 

The reader may wonder whether a solution in the sense of distributions (possibly 
defined by the identity in ^2 Lemma 11-3]) would not be good enough for the 
purposes of this paper. It turns however out that the information yielded by the 
order parameter x in Definition 13.21 carries information that is essential in what 
follows. Incidentally, x may be different from X{ti>o} (see ^ Remark 4.1]). 

Lemma 3.4. Let {u,x) o, solution in the sense of Definition VS. 'A and suppose that 
for some (ccg, to) in the set of definition and for some sequence r,„ 0, m — > oo 

Ur^{y,s) :— ^ ° — ^ — — — ^ locally in {y„ < 0} x (— oo,0) as m oo 
rm 

and 

Xr,„ (y, s) -.^ x{xo + rmy, to + r™^s) -> a. e. in {?/„ > 0} x (-oo, 0) as m ^ oo . 
Then for some 6 > 0, u is caloric in Qs{xo, to) md satisfies 

u^O in Qj{xo,to) ■ 
Proof. The assumptions imply by Definition 13.21 1) that 

a.e. in {x„ > 0} x (— oo, 0) as jti — > oo . 
Moreover, they imply by ^2 Proposition 10.1 2)] that the density 

*(.o,to)(0+)e {0}u{i/„} , 
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where Hn is tlie energy of the half-plane solution defined in Section 10]. In the 
case 

*(.o,t„)(0+) = 

we obtain from Jl] Proposition 10.1 2)] immediately the statement of the lemma. 
In the case 

*(xo,to)(0+) = -ff- 

it follows from [111 Proposition 10.1 1)] that the limit of Ur^{y, s) as m — *■ oo must 
after rotation be the half-plane solution max(— 2:„,0), a contradiction to the limit 
of Ur^ being 0. □ 

4. Flatness Classes 
Definition 4.1. Let < (t+, (t_ < 1 and r > 0. We say that 
u e F{a+, (T_, r) in Qp in direction e„ 

if 

(1) (it, x) is a solution in the sense of Definition 13 . 21 in a domain containing Qp. 
(2) 

{o,P^)Gd{u>o}, 

u{x, t) = x{x, t) = when {x, t) G Qp and Xn > c+p, 
xix, t) = 1 and u{x, t) > — (a;„ -|- cr^p) when {x, t) G Qp and a;„ < —U-P ■ 
(3) 

|Vu| < 1 + T in Qp. 

When the origin is replaced by (xojto) a-nd the flatness direction e„ is replaced 
by u then we define u to belong to the flatness class i^(cr_|_, (t_, r) in Qp{xo,to) in 
direction ;/. 

5. Flatness on the side of {% = 0} implies flatness on the side of 

{X>0} 

The aim of this and the following sections is to draw information from proper- 
ties of an inhomogeneous blow-up limit. One of the central problems when using 
blow-up arguments is ^^not-strong convergence^^ or ^^energy /oss" in the limit. Here 
we avoid those problems by working with uniform convergence (not some Sobolev 
norm). The approach is based on a powerful idea by H.W. Alt-L.A. Caffarclli who 
used "flatness on the side of {u = 0} implies flatness on the side of {u > 0}" to 
prove uniform convergence to an inhomogeneous blow-up limit (cf ^ Section 7]). 
In this section we extend their result to a weaker class of solutions and to the par- 
abolic case, using results in |11) . 

The following Lemma is the parabolic version of Lemma 4.10]. 
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Lemma 5.1. Let (u, x) solution in the sense of Definition W^ in a domain 

containing the closure of a non-empty open ball B = {{y, s) : \{y, s) — (j/o, so)| < c} 
such that -B C {x = 0} and B touches the set {u > 0} at the origin. 
Then 

r u{x,t) 

limsup — — 7 r = 1. 

{u>o}3[x,t)^o pardist((a::, t),B) 

Proof. Let Yfe = {uk, Sfe) G R"+^ be a sequence such that 

n y U{x,t) u{Yk) 

£ = lim sup — J-, r — —r = lim 



{u>Q}3{x,t)~.o pardist((x, t),B) fc^oo pardist(Yfe, S) ' 
Set dk := pardist(yfe, B) and let {xk,tk) = Xk E dB be such that pardist(Yfc, Xk) = 
dk- 

We consider the blow-up sequence 

, , u{dkX + Xk,dlt + tk) , ^, , j2_i , ^ ^ 

Uk{x, t) = — ^ -.Xk{x, t) = x{dkX + Xk,dit + tk). 

We know that, passing to a subsequence if necessary, Uk ~^ uq locally uniformly in 
R"+^ and Xk Xo weakly-* in L5^^(R"+^) as A; ^ oo. Also, after a rotation and 
translation, the scaled B converges to {.t„ > 0} and ((yfc — Xk)/dk, (sfc — tk)/d1) 
(5,r) e dQi{0) as fc ^ oo. The limit function uq satisfies 

Alto - dtuo = in R"+i n {uq > 0} , 
uo(^, t) = £ and 
uo{x,t) = in {xn > 0}. 
By the definition of the limit superior we know also that 

uoixjt) < —ixn in {Xn < 0}. 

The strong maximum principle (applied to uo{x,t) + ixn) tells us therefore that 

uo{x, t) = £max(— x„, 0) for t < t. We have to show that I — 1. 

In the case £ > we obtain from the fact that (m, x) is a solution in the sense of 

Definition 13.21 that xo = 1 in {xn < 0} n < r}. Furthermore, we infer from the 

assumption that xo = in {.t„ > 0}. But then (uo(-, t+r), xo(', ^+''')) is in < 0} a 

solution in the sense of Definition l3 . 21 whose energy M{uq{-, i + r), Xo('5 ^ + ''')) — Hn 

(cf. [Ill Section 10]), whence j^fl Proposition 10.1] implies that (. = \. 

In the case £ = we apply Lemma [3. 41 to obtain for some (5 > that u is caloric in 

Qs and satisfies 

w = in . 

As {u ~ 0} contains _B, u being caloric in Qg and therefore analytic with respect 
to the space variables implies 

w ~ in 

for some 5i > 0. This is a contradiction in view of the origin being a free boundary 
point. □ 
The following theorem extends Lemma 7.2]. 
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Theorem 5.2. There exists a constant C G (0, +oo) depending only on the space di- 
mension n such that ifu € F{a, 1, a) in Qp then u e F(Ca, Ca, a) in Qp/2(0, y-m 0)) 
for some |y„| < Ca. 

Proof. The idea is to touch the boundary d{x = 0} with the graph of a C^-function, 
to apply Lemma Is. II and to proceed then with a Harnack inequality argument. 
Step 1 (Touching d{x = 0} with a smooth surface): 

Rescaling Up{x,t) := "(p^^p ^ Xpi^, t) := we may assume that p = I. 

Let 

_ / eMj^Ww^)^ l^'l' + 1^ - 1|< 1/16, 
0, else 

and let s be the largest constant such that 



77(x',t) 



Qi n {u > 0} C {{x,t) e Qi : Xn < a - sr]{x',t)} =: D. 

This implies that there exists a point {xq, to) Z E dD n d{u > 0} n {i > 15/16}. 
As (0,1) is a free boundary point, we know furthermore that s < a. 
Let us also define the barrier function v by 

Aw — dtv = \ii D , 

w = on dD n Qi and 

v = 2a~Xn ondDnd'Qi. 

Note that this implies that —a < v + Xn < Scr. 

Since |Vu| < 1 + cr we also obtain that v > u on dD and thus, by the maximum 
principle, that v > u in D. As 77 is a C^-function, the assumptions of Lemma f5.il 
are satisfied at Z. Therefore 

(5.1) 1 < limsup— -^^^4w- < -d^Z), 

[x,t)-.z pardist((x,t),B) 

where v is the outward space normal to dD at Z . In order to obtain an estimate 
from above we define 

F{x, t) = 2(7 — Xn — v{x, t). 

F is caloric in D and satisfies < F < a. Since is a regular parabolic domain, 
we know from standard regularity theory for parabolic equations that sup^j |Vi^| < 
Cicr. Therefore 

-dnV{Z) = 1 + dnF{Z) < 1 + CiCr. 

By the flatness assumption we know that v is close to e„. More precisely, 



tfYihl^^SI^\<vro\w,\s. 

V/S2|V7?|2 + 1 I - 



Thus 

-d^v{Z) = -VviZ) ■ {ly - en) - dnv{Z) < 1 + C\a + VTO\Vt]\\Vv{Z)\s < l + Cacr 
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From inequality H5.1|l we infer that 

(5.2) 1 < -duv{Z) < 1 + C2cr. 

Step 2 (Harnack inequality argument): 

As wc know already that v is cr-closc to ^x„, it is sufficient to show that u is cr-close 
to V on the set {{x,t) : Xn = —3/4, \x'\ < f/2, t < 3/4}. Once this is done, we 
may integrate u in the a;„-direction to establish the lemma. 

In order to prove the a-closeness we define for = (7, r), r e (—1, 3/4), |7'| < 1/2 
and 7„ = —3/4 the function o;^ by 

Aiv^ - dtUJ^ = in £> n {t > t} 
io^^-Xn on Bi/g(7) X {t = r} 

= on the remainder of the parabolic boundary of 13 n {t > r}. 

By the Hopf lemma we have 

dyUJ^{Z) < -a < 

uniformly in ^. 

We would like to show that u > v — C^aXn- The trick is to compare u to v ~ Kauj^ 
on the set Bi/^(^) x {t = r} and to use the information on the normal derivative 
of w at Z to prove that if K is large, then u > v — Kulj^ for at least one point in 
^1/8(7) X = ''■}• More precisely: 

Assume that u < v — Kau;^ in Bi/g,{j) x{t = t}. Then u < v — Kaui^ in Dr\{t > r}. 
Consequently, we obtain from inequalities (15.11) and 1)5.2(1 that 

1 < -d^v{Z) + Kad„uj^{Z) < 1 + C2(t - Kaa. 

This yields a contradiction when K is large enough, say K = ^C^jct- Thus u{X() > 
v{X^) — K(juj^{X^) for at least one point G i3i/g(7) x {< = t}. 
On the other hand, v — u>Q. Therefore we can apply the Harnack inequality and 
deduce that 

(i;-u)(|)<C3 inf {v-u)<CAa, 

Qi/8(C+(oa/32)) 

for every |e {(x',-3/4,i) : \x'\ < 1/2, -1 < t < 1/2}. 

This implies that —3/4, t) > 3/4 — C^a in the above region. Integrating in the 
e„ direction and using the assumption |Vu| < 1 + cr yields the estimate 

u > —{xn + Cga) in {—3/4 < Xn < — cr} x Q'1/2 

By our initial assumption we also know that m = in {3/4 > Xn > cr} n Q'lj^- 
Translating (u,x) in the e,i direction so that the point (0,1/4) e d\u > 0} and 
using X ^ X{u>o} of Definition 13.21 1) we obtain the statement of our theorem. □ 
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6. Inhomogeneous Blow-up 

In this section we consider inhomogeneous scahng of the solution and the free 
boundary. The foUowing lemma is our version of ^ Lemma 7.3] 

Lemma 6.1. Suppose that Uk G F{ak,<7kiTk) in Qp^, that cr/j — > and that 
Tk/cr^ ^ 0, and define 

/+ {x', t) := sup{h : lim sup r"^"^ f x > 0}, 

fk{x',t) := M{h : limsupr-"-^ [ X > 0}. 

r-*0 JQr(pkx',akPkh,plt) 

Then, as a subsequence k — > oo, fj^ and fj7 converge in L^^{Q'i) to some function 
f , and f is continuous in Q'^. 

Proof. Rescaling as before we may assume that = 1. Let 

Dk := {{y',h,t) : lim sup r-"-^ f x > 0} . 

r^O JQr-iv'-Xrkh.t) 

We may assume - passing if necessary to a subsequence - that Dk converges with 
respect to the usual (not the parabolic) Hausdorff distance as fc ^ oo. Let us define 

f{x\t):^ limsup 4^(y',s), 

(y'.s) — ^{x\t),k — *oo 

where we take the limit superior with respect to the above subsequence. For every 
{yo,to) there exists then a sequence {y'^^tk) — > iyo,to) such that fkiVk^tk) 
fiy'oito) as k —+ oo. By definition / is upper semi-continuous. Therefore we obtain 
for e > and sufficiently large k that 



(Q'M^tk) X [f+{yltk) + S,^))nDk = (b. 



Consequently Uk e F{ak§, 1, t^) in Q^iyk, c^kft {Vk^ik), tk)- Applying ThcorcmE3 
to Uk we deduce that 

Uk{x,t) > -{xn + CakS/2) for {x,t) G Qe/2iy'k,(^kfk{y'k,tk),tk) ■ 

In terms of f^ and /jT this yields fk{y',t) > fkiv'k^tk) - CS in Q'^/^(y'k,tk)- It 
follows that liiiifc^oo fk{y\ t) — f{y\t); that f^ and f^ converge locally uniformly 
and that / is continuous. □ 
The next Proposition follows the lines of Lemma 5.7]. 

Proposition 6.2. Suppose that the assumptions of Lemma \6. 1\ are satisfied and 
that k is the subsequence of Lemma \6.1\ Then 

( f I, _ Ukjpkx' , Pkh, pjt) + pkh 

(7k 
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is for each d G (0, 1) bounded in Qi-s H {xn < 0} (by a constant depending only on 
S and n) and converges on compact subsets of Qi in to a caloric function w. 
Moreover, w{x' ,h,t) is non- decreasing in the h-variable in and 

lim Wk{y,s) = f{x' ,t) \ 

here f is the function defined in Lemma \6.1\ 

Proof. Rescaling as before we may assume that pk = 1- 

The function is caloric in Qi f] {h < — CTfc}. Using Dcfinition l4.1l 3). we obtain 
that 

Uk < ~x„ + 2ak in Qi n {x„ < 0} , 

implying that Wk < 2. From Thcorcm l5.2l and Definition l4.1l 3) we infer furthermore 
that Uk{x,t) > — (x„ + CgOk) for (a;',x„,i) G Qi-s H {a-„ < 0}, implying that 
Wk > -Cs in Qi^s n {xn < 0}. 

By Definition 14.11 3'! and the assumptions, |Viifc| < 1 + o(cr^.). Consequently, 

(6.1) —dhWk < l^^^^-l 5- < — ^Oasfc^oo. 

In the remainder of the proof we will show that w attains the boundary data / as 
h ^ 0. First, we show that for fixed L G (1, +oo) 

(6.2) Wfc(x', CTfc/i, t) - f^{x',t) uniformly in Q[_s x {-L < h < 0} 

as fc — > oo. An estimate from above can be obtained easily from inequality (|6.1|l : 
Wk{x,hcrk,t) - fk{x\t) < Wk{x,akfk{x,t),t) - fk{x,t) + {fk{x ,t) - h) — 

CTk 

< {1 + L)— ^ as fc oo . 

This establishes an estimate from above. In order to derive an estimate from below 
we use Theorem 15.21 Consider a sequence of points (a;'j,,tfe) G Qi^g and fixed 
S G (4,+oo). Then 

Uk e F{ak, l,Tk) in Qsak{x'k,<Jkfk{xk^tk),tk) 

for 

^k = ^ sup {f+{x',t)~f+{x',,tk)). 

From the uniform convergence of to the continuous function /, we infer that 
CTfc ^ as fc ^ oo. Now by Theorem 15. 21 

Uk G F{Cd-k,Cak,Tk) in Qscrf,/2ix'k,a-kfk{x'k,tk) + CSd-kO/2,tk), 

where ct^. ~ max((Tfe, Tk) and 6 G [0, 1]. 
Thus for h G (max(-L, -5/4), 0) 

Ukixk + haken,tk) > -(Jk {h - fk {x'ki'tk) + CakS) . 
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Consequently 

f , , . . Ukjxk + h(Tken,tk) + hak r+f / . ^- c- 
and lEiSl holds. 

To establish liniQ-g(y^-,_^(^, p j^^g,^ Wfe(y, s) = f{x',t), we need to extend the 

convergence (|6.2ll to larger values of h. To this end, we define the barrier function 
Ze by 

Aze - dtZi, = in Q^_g , 

Ze = ge on 5'(5i-5 n {/l = 0} , 

= inffc infg- Wk on d'Qi-s n {/i < 0} , 

where ge G C°° and f — 2e < < f — s. By Ht).2|l we know that Wk > on 
9'{Qi-s n {/i < — -LfTfc}). From the comparison principle it follows that Wk > z^ in 
Qj^^^- n {/i < — LfTfc}. Thus, by local boundary regularity for solutions of the heat 
equation, liminfQ-_^^3(^_^)_^(^, ^ Wkiy,s) > ge{x' ,t) > f{x',t) - 2e. 

The opposite inequality follows from a similar argument, comparing Wk to the upper 
barrier z defined by 

Az. - dtZe = in Q^_g , 

Ze = ge on d'Qi-s n{h = 0} , 

Ze ~ supj, supg- ^ Wk on 3'(5i-5 n {/i < 0} , 

where 5e e C°° and / + 2e > > / + e. □ 

7. Scaling discrepancy and C°°-regularity of blow-up limits 

In order to obtain "bctter-than-Lipschitz" -regularity of the inhomogencous blow- 
up limit /, H.W. Alt-L.A. Caffarelli used the non-positive mean curvature of d{u > 
0} at singularities. The analogue of the non-positive mean curvature property can 
still be proved in the time-dependent case, however that path leads to problems in 
the sequel. Therefore we replace it by a scaling discrepancy argument which gives 
hope to be applicable in more general situations. We obtain C°° -regularity of /. 

Proposition 7.1. Suppose that the assumptions of Lemma \6. 1\ are satisfied and 
that k is the subsequence of Lemma \6.1\ Then dnW ~ Q on Q\i2 *^ sense of 
distributions. 

Proof. Rescaling as before we may assume that pk = 1. 

In what follows, g{x\ t) = 8(|a;'p-|-|t|)— 4. Note that / > g in Q'if2- ^ct us introduce 
the following notation: Z shall be the set {(a;',a;„,t) : {x' ,t) g Q'i,Xn G R}. Given 
a function <j) : Q'l ^ R, we divide Z into the three parts 

Z+{(t))^{{x,t)^Z : Xn>(j){x' ,t)}, 

Z-{<t>) = {(x,t) eZ: Xn< 

Z\(l)) ^ {{x,t) e Z : xn^^{x\t)}. 
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Moreover let ^ be defined by fi{A) := 'H"^'^{A {s = t}) dt for any Borel set 
A C R"+^. Adding an arbitrarily small constant to the function g, we may assume 
that fi{Z'^{akg) CiRk) = for all k; here Rk is the regular part of the free boundary 
d{uk > 0} introduced in ^] Proposition 9.1], i.e. 

Rk{t) := {x G d{uk{t) > 0} : there is ^^^{x.t) S 9i?i(0) such that Vr{y,s) = 

— '- ^ — > max(— y • i^/fj. (x, t), 0) locally uniformly in {y,s) G R"^^ 

as r — > 0} . 

Last, we define Ek ■= {uk > 0} n Z^{(Jkg) and Sj, ;= {{x' ,t) : [x' ,akg{x' ,t),t) G 
{uk > 0} n Z}. By the choice of g we know that the limit inferior of the sets 
contains Q\i2- 

We will deduce the result from the following three claims. 
Claim 1: 



ti{Z+{akg) n < - / {dnUk + l)dx'dt + + C^ai. 

Claim 2: 

- C2<yl < ti{Z+{akg) n Rk). 

Claim 3: 

/ \dnWk{x' ,akg{x' ,t),t)\ as fc ^ oo . 

Proof of Claim 1: By the representation theorem Lemma 11.3] we know that 
for non- negative (j) G , 

(7.1) [ [ (P dn''-'^ dt < - [ {Vuk-Vcp + dtUkcf)) dxdt . 
Letting (j) Xz+{crkg)XQ2 the inequality (|7.1(l becomes 

(7.2) fi{Z+{akg)nRk)^ [ f dH'^-Ut 

J-oo J Rk{t)l^Z+(rTkg) 



< / Vwfc ■ V dx dt — / dfUk dx dt , 

where u is the outward unit space normal on dZ^ [ukg). Since 

V = ^ =(ak^'g, -1) , 

we obtain 

fi{Z+{akg)nRk) < [ {\/uk){x',akg{x',t),t)- {akVg{x',t),-l)dxdt 

dfUk dx dt . 

{«fc>o}nz+((jfcg) 
Let us rewrite the integral 

{^Uk){x',akg{x',t),t) ■ {akVgix',t),-l)dxdt 
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ak{Vuk){x ,iJkg{x , t),t)-V' g{x ,t)~{dnUk{x , augi^x ,t), t) + l)dx'dt + C{Y.k) 

(7kUk{x\(7kg{x , t),t)A'g{x', t) - ak^d„Uk{x', 0-^.9(2;', t), t)\Vg{x', t)|^ 

- (dnUkix', akg{x', t),t) + 1) dx' dt + C\T.k) 

+ / akUk{x\akgix',t),t)d,jg{x',t) dn"~^ dt , 

where 77 is the outward space normal on Since Uk = on SSfc, the last integral 
is 0. 

Moreover, A'g = 16 and Uk < C^dk on {x' , g{x' ,t),t), implying that 
(VMfc)(x',afc(7(x',t),i) • {akV'gix',t),^l)dxdt 

(dnUkix', cTkg{x', t),t) + l)dx'dt + £"(S]fc) + dal . 
By the definition of Wk this tells us also that 

(7.3) / {Vwk)ix\ akg{x', t), t) ■ {ak^' gix\t), 0) dx dt ^ ask ^00 , 



a fact that will be used later on. 

Last, integration by parts of the last term in (|7.2() with respect to the time variable 
yields 

- / dtUk dx dt < CsCTfc^ . 

J{uk>0}nZ+{<jkg) 

Combining the above estimates we obtain Claim 1. 

Proof of claim 2: With the outward space normal on the boundary of Z^{akg) 

1 



V. 



and with the outward space normal h'ji^, on the regular boundary of Ek we compute 



(7.4) fi{Z+{akg)nRk)> 
f-i 



Vg^ ■ vr^ dH" ^ dt 



1 J Z+{<jkg)nRk{t) 
1 

div i/g, dH"~^ dt 

1 JEknZ+{at:g) 

1 

I/g,, • Vg^dW'^ dt. 



The normal satisfies 



-1 JdZ+{akg)nEt, 



div i>„, > — , > — CfiCrfc. 

Inserting this estimate for the divergence into H7.4(l yields 
(7.5) KZ^i^kg) n Rk) > ii{dZ+[akg) n Ek) 
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CeCTfc dW'-^ dt 

1-1 JEkr\Z+(okg) 

> fi{dZ+{akg)nEk) - Oral ; 
the last inequality follows from the fact that the width of the set Ek is of order 
0{ak)- As the area of dZ^{akg) n Ek) is greater than that of S^, the statement of 
Claim 2 holds. 

Proof of Claim 3: From Claim 1 and Claim 2 we infer that 

-CsCTfe < - / idnUkix',cTkg{x',t),t) + l)dx' : 

But since G F{(Tk, (Jk,Tk) and Tfc/cr^ — > as fc — ^ oo, it follows that 
a„ufc + l > -iVwfel + 1 > -o{al). 

Consequently 



dt 



\dnWk{x' ,(Jkg{x\t),t)\ 



dnUk{x' ,(Tkg{x' ,t),t) + 1 



0-fc 



< I 2 max 



dnUk{x' ,akg{x' ,t),t) + 1 



,0 



dnUk{x' ,(7kg{x' ,t),t) + 1 



< C<j<Tk as fc — !■ cx) , 



and Claim 3 is proved. 

Proof of the Proposition: Let C G Cq(Qi/2). From Claim 3, from the fact that Wk 
is caloric in Z~{akg), from H7.3|l and from a standard energy estimate for caloric 
functions we infer now that 



0(1) 



C dnWkix, (Jkg[x\ t), t)Vn 



{dnC d„Wk - C^'wk + CdtWk) 



= 0(1) + / {^r^C dnWk + V'C-V'u.fc - WkdtO 

—* I {dnC dnW + V'C • V'w — w dtQ as fc ^ oo ; 

here v is the outward unit space normal on dZ^{akg)- It follows that dnW = on 
Qi/2 sense of distributions. □ 

Corollary 7.2. Suppose that the assumptions of Lemma \6.1\ are satisfied and that 
k is the subsequence of Lemma \6.1\ Then f € C'°°((5i/2)/ moreover, 

Qa+k j 



<C{n,\a\,k) 



dx°'dt^ 

in Qi/i for any fc G N and multi-index a € N". 

Proof. Since dnW = on Q'j^^j ^^"^ sense of distributions we may reflect w to a 
caloric function in Qi/2- As / = '^\q'i ^'^'^ ll'"^IU°°(Q3/4) — ^W) (see Proposition 
I6.2|l . the result follows from standard regularity theory of caloric functions. □ 
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8. Flatness improvement and regularity 

Concluding regularity is then a standard procedure. The following Lemma l8. II 
Lemma [8.31 and Theorem 18.41 extend Lemma 7.9, Lemma 7.10 and Theorem 8.1 in 
Finally, we apply Theorem 18.41 to regular free boundary points, i.e. points in 
the set R defined in m Proposition 9.1] (or the proof of Proposition l7.1|l to obtain 
that R is open relative to d{u > 0}. 

Lemma 8.1. Let 9 G (0, 1). Then there exists a constant ae > depending only 
on and the dimension n such that if a < ae, t < aea"^ and u G i<"((T, ct, t) in Qp 
in direction rj, then 

ueF{ea,l,T) in Q^(^n)ep{^V,0) 

in direction rj for some i? G [— cr, cr] and some rj satisfying \rj ^ ri\ < C{n)a. Here 
c{n) > and C{n) < +oo are constants depending only on the dimension n. 

Proof. We may rotate the coordinate system so 77 = e„, and we may assume that 
/3 = 1. By a contradiction argument, it is sufhcient to prove the statement of the 
lemma for Wfc as in Lemma l6. II and every large k. 
First, observe that by Corollarv l7.2l 

f{x',t) < f{o,o) + e- x' + c{\x'\^ + \t\) in q;/4, 

where £ is the space gradient of /, \i\ < C and C depends only on the dimension 
n. Thus 

fix', t) < /(O, 0)+i.x' + ^^m Qg^^^c) ■ 
It follows that for large k the function in Lemma [6.11 satisfies 

ftix', t) < /(O, 0) + £.x' + e^ in ge/(4C) • 
This means that Uk G F{9a, l,r) in (56(/(4C)(0, /(O, 0), 0) in the direction fj, where 

^/^TK^■ 

The lemma follows. □ 
Lemma 8.2. Let u he a solution in the sense of Definition \3.'^ Then 

max(|Vw|^ - l,0)(a-,i) ^ as < pardist((2;, t), {li = 0}) ^ . 
Proof. Consider a sequence {u > 0} 9 {xk,tk) — > (xo,to) such that 

l<e:= limsup \Vu{x,t)f = lim \Vu{xk,tk)\'^ . 

Setting rfc := pardist((xfc, t^), {w = 0}), the blow-up sequence 

u{xk + rky,tk^rk^s) , , 2n 

Uk{y,s) := ,Xk[y,s) := x{xk +rky,tk + rt s) 
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converges to a solution (uo,Xo) in the sense of Definition 13.21 satisfying uq > in 
Qi, |VmoP < £ and |Vuo(0)p = £. The strong maximum principle implies that 
UQ{y, s) = £max{y • e, 0) in {y • e > 0} H {s < 0} for some e e dBi. 
From 11, Theorem 11.1] we infer that 

{y • e = 0} n {s < 0} C S** 

up to a set of vanishing -measure, where 

E**(t) {x e d{uo{t) > 0} : there is e{x,t) e (0,1] and ^(x,t) e dBi{0) such 

that — — 9{x,t)\y ■ ^{x,t)\ locally uniformly 

r 

in {y, s) e R"+i as r ^ 0} . 
However 6{x, t) G (0, 1] contradicts £> 1. □ 

Lemma 8.3. For every 6 G (0, 1) there exist ug > Q and cg G (0, 1/2) depending 
only on 6 and the dimension n such that if u G _F((t, 1,t) in Qp in direction rj with 
cr < ag and r < aga^ then u G F{9a,0a,6^T) in Qcgpiy,0) in the direction rj for 
some y,rj satisfying \rj ~ ri\ < C{n)a and \y\ < C{n)a. Here C{n) depends only on 
the dimension n. 

Proof. We may assume that p = 1. 

From Lemma 15.21 we infer that u G F(Ca,C(j,T) in Qi/2(2/,0) in direction i] for 
some y G Bca- Consequently we may apply Lemma |8. II to deduce that for some 
9i to be determined later, u G F{C6i<7, 1, r) in Qc(n)ei iVi 0) in the direction 77 such 
that \rj — fi\ < Ca and |y — yl < (C + 1)(t < 1/2, provided that ag has been chosen 
small enough in terms of 9i . 

In order to be able to continue we need to show improvement with respect to the 
r-variable. To that end, observe that U = max(|Vu| — 1,0) is by Lemma 18.21 a 
continuous subcaloric function in Qi with boundary values less than tx{u>o} ^ 
TX{xn<cr}- We may therefore compare U to the caloric function with boundary val- 
ues TX{x„<a}- It follows that < U < (1 — ci)r in Q1/2 for some ci > depending 
only on the dimension n. Thus u G F{C6ia,l, {1 — ci)t) in Qc(n)ei(y5 0) ij^i 
direction i]. Choosing 9q := ^/l — ci and 9i := 9o/C we obtain u G F{9Qa, 1,9qt) 
in Qc^Ooiy^ 0) in the direction fj such that |?7 — ?7| < Ccr, where C2 G (0, 1) depends 
only on the dimension n. 
Iterating this process we see that 

u G F{9'^a, l,6'g"V) in ^(c^eo)- (y™, 0) in the direction f]m 

where \r] - f]ra\ < C{n)aYJ^=o and \ym\ < C[n)aY^jLo{c20oy ■ 

Applying once more Lemma and choosing 6*0 := 9^ /C we obtain the statement 

of the lemma. □ 
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Theorem 8.4. There exists a constant do > such that if u ^ -F((t, 1,t) in 
Qp{to,xo), cr < ctq and t < aocr^ , then the topological free boundary d{u > 0} 
is in Qp/4{tQ,XQ) the graph of a C^^"-"^ -function; in particular the space normal is 
Holder continuous in Qp/4{tQ,xo). 

Proof. Using Lemma l8 . 31 inductively we see that 
(8.1) 

u e F{0^a, 6^ a, O'^^t) in Q^^ kp{y, s) in the direction 77*^ 

where \t - v\ < Cin)aj:fz',m' and \t ~ y\ < C(n)aE ■ro(2ce/20)^' , 

provided that (y, s) G Q 1/2(^0, 2^0) H d{u > 0}, 9 < 1/4 and 

fxo <min(l/(4CH),a(,/2/2) ; 

here we sacrificed some flatness in order to keep the original free boundary point 
(y, s). We obtain existence of the outward space normal v on Qi/2{to,xo). More- 
over, 1/ satisfies by (j8.1(l 

which implies Holder-continuity of □ 

Corollary 8.5. For each point (xo,to) of the set R, the topological free boundary 
d{u > 0} is in an open neighborhood of {xq, to) the graph of a C^~^°''°' -function; in 
particular, the space normal is Holder continuous in an open space-time neighbor- 
hood of (xo, to)- 

Proof. The Corollary follows from |11[ Proposition 9.1] and Lemma [8.21 □ 
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